Global textures and the Doppler peaks 
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Abstract. We review recent work aimed at showing how global topolog- 
ical defects influence the shape of the angular power spectrum of the cos- 
mic microwave background radiation on small scales. While Sachs- Wolfe 
fluctuations give the dominant contribution on angular scales larger than 
about a few degrees, on intermediate scales, 0.1° < <2° , the main role 
is played by coherent oscillations in the baryon radiation plasma before 
recombination. In standard cosmological models these oscillations lead 
to the 'Doppler peaks' in the angular power spectrum. Inflation-based 
cold dark matter models predict the location of the first peak to be at 
£ ~ 220 / "v/Ho , with a height which is a few times the level of anisotropies 
at large scales. Here we focus on perturbations induced by global tex- 
tures. We find that the height of the first peak is reduced and is shifted 
to ^ ~ 350. 



1. Introduction 

The importance of the observations of anisotropies on intermediate and small 
angular scales cannot be over-emphasized. Presently a number of experiments 
scrutinize different regions of the sky trying to uncover the real amplitude and 
structure of the relic radiation. The CMB anisotropies are a powerful tool to 
discriminate among the different models of structure formation by purely linear 
analysis. These anisotropies are parameterized in terms of C^'s, defined as the 
coefficients in the expansion of the angular correlation function 



For Harrison-Zel'dovich spectra of perturbations £{£ + l)Ci is constant on large 
angular scales, say £ < 50. Both inflation and topological defect models lead 
to approximately scale invariant spectra on large scales. The next important 
step towards discriminating between competing models is the measurement of 
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degree scale anisotropies. Disregarding Silk damping, gauge invariant linear 
perturbation analysis leads to (Durrer 1994) 
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where ^ and ^ are quantities describing the perturbations in the geometry and 
V denotes the peculiar velocity of the radiation fluid with respect to the overall 
Friedmann expansion. In Eq.(|2|), D^p = 5pr/ Pr specifies the intrinsic density 
fluctuation in the radiation fluid; below we use Dr = {Spr + Spb)/{pr + Pb), about 

5% smaller than Dg'^ for scales inside the horizon, which are the ones we are 
interested in. 

The Sachs-Wolfe (SW) contribution in the above equation has been calcu- 
lated for both inflation and defect models, yielding mainly a normalization of 
the spectra. We present below a computation for the Doppler contribution from 
global topological defects; in particular we perform our analysis for vrs-defects, 
textures (Turok 1989), in a universe dominated by cold dark matter (CDM). We 
believe that the main ideas of our analysis remain valid for all global defects. 
The Doppler contribution to the CMB anisotropies is approximately given by|^ 
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In the previous formula 7 denotes a direction in the sky 
I time, with r/o and T]rec the present and recombination 



where Xrec = x + 7770. In the previous formula 
and Tl is the conformal time, with TjQ and fjrec 
times, respectively. 



2. Linear theory povi^er spectra 

We now study a two-fluid system: baryons plus radiation, which prior to recom- 
bination are tightly coupled, and CDM. We start off by considering the evolution 
equations in each fluid component as given by Kodama & Sasaki (1984). The 
evolution of the perturbation variables in a flat background, = 1, can be 
expressed as 
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D',-3w<^Dr = {I + W)[3^^ - 3<^' - kVr - I {^^) k-Hl + ^)Vr] 

D', = 3^M'-3<I>'-A;y,-|(^)'fc-i(l + ^)K, 



(4) 



where subscripts r and c denote the baryon-radiation plasma and CDM, respec- 
tively; D, V are density and velocity perturbations; w = Pr/Pr, = v'rl P'r and 

P = Pr+ Pc- 



^We are neglecting here the SW effect, which decays on subhorizon scales like . We also 
neglect the contribution of the neutrino fluctuations. These approximations lead to an error 
of less than about 30% in the amplitude of the first Doppler peak; see discussion in (Durrer, 
Gangui & Sakellariadou 1995). 
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Throughout this section we will normalize the scale factor such that a{rif,g) = 
1, where 7?eg is the conformal time at which matter and radiation densities co- 
incide. Setting r = (SnGpeq/Q)'^^'^ , where peq is the total energy density at 
equality, we obtain a{rj) = (r7/r)(l + jirj/r)). It is clear from this expression 
that for early (late) conformal times we recover the radiation (matter) dominated 
phase of expansion of the universe. 

Since the energy density of baryons and CDM go like oc a~^, the ratio Pb/Pc 
is const and equal to Qb- Therefore {pb/ Prad) — f^^a. This leads to 

wc^^-il + ^Ba)-^ ; cl^'^-{l + -^^Ba)-\ (5) 

The only place where the seeds enter the system (Q) is through the potentials 
^ and These potentials may be split into a part coming from standard 
matter and radiation, and a part due to the seeds, e.g., ^ = ^(c,r) + ^seed where 
^ seed-, ^seed are determined by the energy momentum tensor of the seed; the 
global texture in our case. Having said this, one may easily see how the seed 
source terms arise. Long but straightforward algebra leads to two second order 
equations for Dr and Dc, namely 

a' a' 

D'; + -[l + 3cl-6w + F-^pM PcF'^il + w)D'^ 

a a 

+ATiGa^[pr{3w'^ -Sw + Qcl - 1) - 2F"^wpc(Pr + Pc) 

1.2 

+Pc(,9cl-7w) + ^-^ci]Dr-4TTGa^Pc{l + w)D, = (l + w;)5; 

D'c +-[! + (! + w)F-^pr{l + 3c1)]D'^ - -(1 + 2>cl)F-^ prD'^ 
a a 

-^■KGa^PcDc- ^'KGa^Pr{l + 'icl)[l-2{pr + Pc)F-^w]Dr = S, (6) 

where F = /c^(127rGa^)~^ + pr{l + w) + pc and S denotes a source term, which 
in general is given hy S = A-KGa^{p + ?>pY^'^'^. In our case, where the seed is 
described by a global scalar field we have S = 8iTG{(j)')'^ . From numerical 
simulations one finds that the average of over a shell of radius k, can be 
modeled by (Durrer & Zhou 1995) 

with fj denoting the symmetry breaking scale of the phase transition leading 
to texture formation. The parameters in are A ~ 3.3, a ~ — 0.7/(27r) and 
P ~ 0.7/(27r)^. On super-horizon scales, where the source term is important, 
this fit is accurate to about 10%. On small scales the accuracy reduces to a 
factor of 2. By using this fit in the calculation of Dr and Dc from Eqs.(|6|) we 
eff'ectively neglect the time evolution of phases of {4>')'^', the incoherent evolution 
of these phases may smear out subsequent Doppler peaks (Albrech et al. 1995), 
but will not affect substantially the height of the first peak. 
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3. Choosing initial conditions 

In defect models of large scale structure formation one usually assumes that 
the universe begins in a hot, homogeneous state. All perturbation variables are 
zero, and the primordial fields (e.g., the global field (/> in our case) are in thermal 
equilibrium in the unbroken symmetry phase. When the phase transition occurs 
the original value (p = becomes unstable and (p assumes a value which lies in the 
sphere of degenerate minima determined by < (/>^ >= r/^. Causality requires the 
correlation length to be bounded by the size of the horizon. On super-horizon 
scales the correlations vanish and therefore the defect stress energy tensor has a 
white noise power spectrum for hrj « 1. After the phase transition, the defect 
energy rapidly enters scaling, p^^g/ oc r]~'^. Adopting this white noise behaviour 
of the field on super-horizon scales, one can simulate (p on scales much larger 
than the horizon scale at the phase transition. The only condition is that the 
scales considered are super-horizon at the time one starts the simulation. This 
is very important, since the scales relevant for large scale structure are more 
than 50 orders of magnitude larger than the inverse symmetry breaking scale 

Going back now to our set of Eqs.(P) we will specify our initial conditions 
as follows: for a given scale k we choose the initial time r/j„ such that the 
perturbation is super-horizon and the universe is radiation dominated. In this 
limit the evolution equations reduce to 

,,2 A Ae ,,3,3,3 Ae 

D'--Dr = -— ; D' +-D',--Dl--^Dr = — , 8 

with particular solutions Dj. = —^eArj^^'^ and Dc = — |eAr/^/^. In the above 
equations we have introduced e = AirGfj'^, the only free parameter in the model. 
We consider perturbations seeded by the texture field, and therefore it is incor- 
rect to add a homogeneous growing mode to the above solutions. With these 
initial conditions, Eqs.(^) are easily integrated numerically, leading to the spec- 
tra for Dr{k,r]rec) and D'^{k,r]rec)- (see Durrer, Gangui & Sakellariadou 1995 
for details). 



4. Angular power spectrum from global textures 

In this section we calculate the expression for Ci as a function of the baryon- 
radiation power spectrum and its derivative, evaluated at recombination, Tfy^^. 
We express the Fourier transform of the velocity perturbation in the baryon- 
radiation component as V{k) ~ —ikD'^{k)/\k'^{l + w)]. We then Fourier trans- 
form the Doppler contribution to the anisotropics given by Eq. (|3|) and perform a 
standard analysis to get the Doppler contribution to the angular power spectrum 
Ci. This is given as a function of the spectra for Dr{k,r]rec) and D'^{k,r]rec) as 
follows 

^Vr(A:,t,ec)Pi|(A:to) + (1 + t«)"'|^;(^,trec)|'(j;(fcto))' 



Ci = - I dk 
vr 
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• (9) 



As we mentioned earlier, the source term in Eqs.(|6|) was estimated from nu- 
merical simulations of the evolution of the defect field, which give a fit for the 
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Doppler peaks from global textures 




multipole number / 



Figure 1. The angular power spectrum for the Doppler contribution 
to the CMB anisotropics is shown in units of (thick solid curve). 
Contributions from the baryon-radiation density perturbation (dashed 
curve) and its derivative (dotted curve) from Eq.(9) are also shown. 

absolute value of but do not specify the phases. The crossed term missing in 
Eq.(^), cannot be determined with this method. However, due to phase aver- 
aging, we expect it to be diminished substantially and therefore we neglect it. 
Integrating Eq.(^), we obtain the Doppler contribution to the CMB anisotropics 
shown in Figure 1. 

5. Discussion 

As we see from the figure, the angular power spectrum + l)Ci yields the 
Doppler peaks. For I < 1000, we find three peaks located at £ = 365, £ = 720 and 
£ = 950. Silk damping, which we have not taken into account here, will decrease 
the relative amplitude of the third peak with respect to the second one; however 
it will not affect substantially the height of the first peak. The integrated SW 
effect, which also has been neglected, will shift the position of the first peak to 
somewhat larger scales, lowering £peak by (5 - 10)% and possibly increasing its 
amplitude slightly (by less than 30%) (Durrer, Gangui & Sakellariadou 1995). 

Our second important result regards the amplitude of the first Doppler 
peak, for which we find £{£ + l)Ci j^^ggg = 5e^ . It is interesting to notice that 
the position of the first peak is displaced by Ai ~ 150 towards smaller angular 
scales than in standard infiationary models (Steinhardt 1995). This is due to 
the difference in the growth of super-horizon perturbations (Sugiyama 1996), 
which is Dr oc r/^/^ in our case, and Dr oc 77^ for inflationary models. 

Let us now compare our value for the Doppler peak with the level of the 
SW plateau. Unfortunately the numerical value for the SW amplitude is uncer- 
tain within a factor of about 2, which leads to a factor 4 uncertainty in the SW 
contribution to the power spectrum: the results are i{£ + l)C£ | ~ 2e^ (Ben- 
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nett & Rhie 1993; Pen et al. 1994) and i{i + 1)Q | ~ Se^ (Durrer & Zhou 
1995). According to the first two groups, the Doppler peak is a factor of ~ 3.4 
times higher than the SW plateau, whereas it is only about 1.5 times higher if 
the second result is assumed. (We allow for about 30% of the SW amplitude to 
be added in phase to the Doppler amplitude of ~ 2.24e.) Improved numerical 
simulations or analytical approximations are needed to resolve this discrepancy. 
However, it is apparent that the Doppler contribution from textures is somewhat 
smaller than for generic inflationary models. 

We believe that our results, about the position and amplitude of the first 
Doppler peak, are basically valid for all global defects. This depends crucially 
on the 1/y^ behavior of ((/>')^ on large scales (cf. Eq. (|^)). The displacement of 
the peaks towards smaller scales is reminiscent of open models, where the first 
Doppler peak is located at Ipeak — 220/^/^Q. A value — 0.3 would produce 
ipeak — 365 for the first peak, as we actually find. However, the amplitude of 
the peaks within open inflationary models is generically larger than the one we 
find for global defects. 

The problem treated in this paper has also been studied by (Crittenden & 
Turok 1995). They obtained the same position of the first Doppler peak but 
with a somewhat higher amplitude. 
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